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PREDICTION  OF  DURABILITY  OF  ELEMENTS  AND  AUTOMATIC 
SYSTEMS  WITH  VECTORIAL 
DETERMINING  PARAMETERS 

-USSR- 

f Following  is  the  translation  of  an  article 
by  G»V.  Druzhinin  (Moscow)  in  law  estiva  Aim, 
deiaiit  N&uk  3SSH,  Otdelenlve  Tekhniehosklkh 
Hank,  Bner.got.lka  1  Avtonatika  (News  of  the 
Academy  of  Sciences  (SSSR,  Technical  Soienoea 
Section,  Power  Production  and  Automation), 

No  S,  Moscow,  1SG1,  pages  1.6 5-170.] 


Durability  means  the  ability  of  elements  and  systems 
to  remain  in  perfect  order  during  storage  or  operation. 
Statistical  studies  of  durability  deal  with  random  change 
processes  in  determining  element  parameters  as  a  result  of 
wear  or  deterioration  (references  1,3),  The  cited  papers 
assume  that  each  element ' Is  characterized  by  a  single 
determining  parameter*  There  are  oases,  however,  when  a 
system  or  element  is  characterised  by-  several  determining 
parameters  <a  woo t or i&l ’.determining  parameter).  It  thus 
becomes  necessary  to  deal,  with  vectorial  random  processes 
of  change  In  the  determining  parameters  and  statistical 
bounds  Of  elements*  The  purpose. of  the  present  article 
is  to  predict  systems  and  element  durability  on. the  basis 
of  the  etate  of  elements  at  the  present  time. 

*-•  Characte  rig  tics,  of  vectorial  deterainin^ 

BKfeflCTL.  a Cjg_boun.de  of  elements.  Any  study  of 
durability  is  almost  always  connected  with  the  extrapola¬ 
tion  of  determining  parameters  and  statistical  bounds. 

At  the  same  time,  the  possibility  of  measuring  determining 
parameter  and  statistical  bound  values  is  extremely  limited. 
This  makes  it  necessary  to  employ  approximate  methods  of 


3. 
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predicting  element  and  systems  durability.  It  is  very 
convenient  to  approximate  vectorial  determining  parameters 
by  means  of  linear  vectorial  random  time  functions  with 
components  of  the  form 

Hv (t)  —  Av  +  Bvt  (v  =- 1, 2, v  (i.j) 

Here  Ay  ie  the  random  initial  value  of  the  \>-th 
component  H  *y  (t)  of  the  vectorial  random  function  H(t); 

Bv  is  the  random  mean  (with  respect  to  time)  rate  of 
change  of  the  \Ath  component  of  H(t>» 

The  formulas  giving  the  mathematical!  expectation*, 
dispersions,  and  correlative  functions  of  the  vectorial 
linear  random  function  HCt)  components  are  analogous  to 
the  corresponding  formulas  for  the  scalar  linear  random 
function  (reference  1).  The  mutual  correlative  function 
of  the  vectorial  linear  random  function  has  the  form 


(tf  O  —  ^0,2,1  +  (£  +  t')^av !>p.  +  Wkbvby,  (i-2) 

Hotting  V  a  in  formula  (1,2),  me  obtain  the 
formula  for  the  correlative  function  of  the  ^-th  component 
of  HCt),  vhile  for  9  -  ^  and  t  »■  t#  we  have  the  formula 
for  the  dispersion  of  this  component. 

Making  use  of  linear  veetorial  random  functions, 

It  is  possible  to  characterise  the  veetorial  determining 
parameters  with  a  relatively  email  number  of  constant 
numerical  characteristics;  this  simplifies  considerably 
the  study  of  element  and  systems  durability.  In  addition 
to  this,  the  use  of  the  indicated  functions  makes  it  possible 
to  extrapolate  the  random  change  proeeeses  in  the  determin¬ 
ing  parameters  and  statistical  bounds  of  elements  on  the 
basis  of  the  minimum  number  of  o rose -sect ions  for  these 
processes  (two).  In  calculating  the  characteristics  of 
linear  random  functions  on  the  basis  of  the  numerical 
characteristics  of  two  cross-sections  of  these  functions, 
each  realisation  is  replaced  with  a  rectilinear  secant; 
there  is  a  statistical  linearisation  of  the  actual  random 


process. 

The  normally  distributed  vectorial  linear  random 
function  is  fully  characterised  by  the  following  constant 
numerical  characteristics  $  the  mathematical  expectations 
for  the  initial  values  of  the  components  m  ,  aa  , ; 
the  mathematical  expectations  of  the  rates  1  ^  of  change 


of  the  components 


®bi*  ®b~,***, 


and  the  central  second-order 


2 


moments 


lot  reduction  of  aiditioaal  a&supptlons  makes 
possible  a  considerable  reduction  in  the  number  of  numerical 
efearaot eristics  describing'  the  linear  vectorial  random  func¬ 
tion  HCt )  *  .  .  . 

Doing  theorera.fi  ©a  the  numerical  characteristics  of 

random  quantities*  it  is  possible  to  derive  formulas  for 
finding  the  numeric®!  characteristics  of  a  linear  vectorial 
random  function  on  the  basis  ©f  the  numerical  characteristics 
of  two  cross-sections  of  this  function..  The  mathematical 
expectations  of  the  initial  values  and  -rates  of  change  of 
vectorial  linear  random  function  H(t)  components  are  ex¬ 
pressed!  by  the  formulas 


M+l" 


m 


m,„ 


m* 


H-r 


■  t. 


< 


(1*3) 

(1*4) 


la.  formulas  (1*3)  and  (1*4}  s.,^  to  the  mathematical 
expectation  of  the  v-th  component  <rO  »•  l»2s,««,h)  of  the 
vectorial  linear  random  function  K(t )  in  the  l«-th  cross- 
see  t  ion* 

The  coupling  moments  of  the  initial  values  . 

f"  V,*l  H-WV«  +  ^VdW  + 

••)•  (v,  pa*  1,2, ..  *,t)  (1*5) 

In  formula  Cl. 3}  3-85  t&®  coupling 

moment  of  cross-section  t±+t  the  v«th  component  of -.the 
linear  random  function  H(t)  and  the  cross-section  of 
the  u—th  component  ©f  the  same  random  function*  The  rest 
of  the  definitions  are  read  off  accordingly*  If  in  formula 
Cl«5)  we  set  V »  a »  we  have  a.  formula  for  the  .dispersion 
of  the  initial  value  of  the  -tb  coaponsat  of  the  linear 

vectorial  random  function  HCt), 

Tha  coupling  moments  of  the  rates  of  change  of  the 

vectorial  linear  random  function  components  are 


t  * 

,_•■•  .  ■  (v,  ji=>  1. 2, . . -tfc) 


'  "t"  "’MV*? 


Tor  v>  we  obtain  a  formula  for  the  diapers  ion 
<fb*  *  fcwW  o£  tke  ^~th  component. 

The  coupling  moment s  of  the  initial  values  of  the 
T  -th  cokponenta  of  the  linear  rectorial  random  function 
H(t)  with  the  rates  of  change  of  the  ^-th  components  of 
this  function  can  be  Aeterainod  f rea  the  formula' 


k  (‘i+r 


Tji  [*i+i^'iviV<‘+x)  ~  f*A*»(t+i)V*+»  ~ 

+  (v>  a  **  i*  2*  •  •  *•  *) 


a.  The  calculation  of  the  durability  _of  elements 
with  rectorial  determining  parameters.  At  a  moment  of  time 
t.  an  elenent  ie  charaoterimed  by  a  system  of  random  quan¬ 
tities  H,^.  let  no  denote  the  mutually  in¬ 

dependent  critical  parameter  values,  at  which  the  element 
Is  considered  out  of  order  (permissible  limits)  by 
uix,  4>h*  Then  the  probability  Q%  that  the  element 

will  be  in  working  order  at  a  moment  of  time  t*  ie  equal 
to  the  probability  of  the  simultaneous  fulfillment  of  the 
inequalities  HXi  >  Hai>  *•••» 

00  00  00 

G*  “J  ^  ^ .  .TJ/.}rfqidri*; .  .drjfc  (2.1) 

where  *i<^\  1*  A  »•  *  *  *  *'<\  >  **  the  probability  density  of  the 
random  rector  with  components  , •«« 

In  practice  it  is  difficult  to  expect  that  the 
rectorial  determining  parameter  of  an  element  will  have  more 
than  three  components*  This  fact  is  engendered  by  mathema¬ 
tical  difficulties  and  the  fact  that  elements  are  usually 
quite  adequately  charaoterimed  by  a  rectorial  random  func¬ 
tion  with  a  small  number  of  component*. 

in  a  number  of  practical  problems  the  components  of  . 
a  normally  distributed  rectorial  random  process  H(t)  are  so 
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w®aklv  correlated  that  they  ca»  be  considered  uncoupled 
&?sd  consequently  independent,  I’hen  integral  <2*l>  disinteg¬ 
rates  into  the  integral  product 

co  co  «o 

Gi  =  6'xG'si .  .  .  Gfe  -  *J/u  (ill)  «?ni  5  hi  (%><*%  •' •••  \  fkt  (n-i)^ih  (2.2) 

to,  »i  “A 

.  \  f  (V\.)  are  the  normal  distribution 

where  *xi^  'I  1^  *  *  *  *  ,xhlv  (hr 

laws  for  the  components  of  the  random  vector  Hi 


i  r 

'<<*»--rrsr“  n — ^rJ 


In  formula  C2*3>»  ®n  is  the  mean  value  of  the 
$-  -th  determining  parameter  of  an  element  at  moment  of 
time  tj  ;  i*  the  mean-square  deviation  of  the  v«th 

determining  parameter  at  moment  of  time  tj.* 

With  normal  distribution  laws  fy^dju')  *  formula  <2*  2) 
can  be  written  in  the  fora 

G,  -  [}•  -  («u)]  *  ■  ■  [i  ■ — ® {Ukl}]  (2A) 


where 


while 


,  1M  r  x»,  1  , 

“PL— rj* 


1®  a  Gaussian  function, 


The  probability  density  for  the  duration  that  the 
element  is  in  working  order  j!u)<t5  is  the  rate  of  change  of 
probability  1-C^  that  the  element  is  out  of  order 


,  ,rt,  _ ;__i  '  (2-»)  . 

/»WN  —  rf«  i*»q 

For  a  vectorial  determining  parameter  with  indepen¬ 
dent  components 

U  (0  Jr-i,  -  €, . . .  Gh  +  -~~t  Gfi-t  ,..Gk  +  <7^  . . .  6h]k  (2.7) 

Since  it  is  usually  the  case  that,  for  elements  with 
rather  short  terms  of  service  the  values  of  *  *  *  *  Gh 


a 


differ  little  from  unity, 


/wU— <2*> 

vwl  1 

Using  formula  (2.8)  it  is  possible.  In  principle, 
to  obtain  an  analytic  expression  for  the  distribution  law 
for  the  time  that  the  element  is  in  good  working  order. 

In  the  general  case,  however,  this  expression  would  be 
extremely  unwieldy.  In  practice,  the  density  of  the  proba¬ 
bility  of  the  time  that  the  element  is  in  good  working 
order  is  conveniently  calculated  by  means  of  the  mean  value 
f  on  each  of  the  intervals  (tif  t1+l>.  In  accordance  with 

(2*  8 )  . . . . ■■■■ . 


Tot  the  normally  distributed  vectorial  determining 
parameter,  formula  (2,9)  is  conveniently  expressed  by  means 
of  the  Gaussian  function 


(2.10) 


viutrs  m is  determined  in  accordance  with  (2.S). 

Thus,  the  approximate  prediction  on  the  basis  of  a 
minimum  number  of  element  parameter  measurements  of  the 
distribution  law  f(t)  for. the  time  that  the  element. ig. In 
working  order  -for  an  element  with  a  normally  distributed 
veotorial  determining  vector  having  independent  components, 
it  is  necessary: 

a)  tc  use  measurements  Of  the  realizations  of  the 
vectorial  determining  parameter  components  of  similar  ele¬ 
ments  at  two  moments  of  time  t  and  t^^  to  calculate  the 

statistical  numerical  characteristics!  mean  values,  disper¬ 
sions,  and  correlative  cross— section  moments} 

b)  using  the  resultant  data  to  find  by  means  of 
formulas  (1.3)-C1.7)  the  characteristics  s»n  Ct)  and 
of  the  vectorial  linear  random  function  H(t>* 

c)  to  make  use  of  formulas  (2.5)  and  (2.10)  to  find 
the  values  of  f^  for  various  time  intervals  (t^,  ti+1>» 

The  calculated  values  of  are  employed  to  construct  a 


6 


hi«to<rru»  which  is  smoothed  out  by  a  continuous  curve. 

3)  The  -.laulatlnu  -  -"-btfltv  of  »*».*«■& 

_..e  .serial  rtwMLMi  ays t emg_ a«t eralnlnK  paraae.oja.^ 

r“««.rsgr=isr»  to  the  i,robie!  a::r;^a 

use  in  various  engineering  applications  can  bn  a_ 

oniv  in  the  ease  of  independent  vectorial  determining  para 
"ter  opponents  fbr  elonents  and  systeus  «  is  then  ponsi- 
ble  to  use  ideas  suggested  in  (reference  3)  for  .Onlar  de- 

termining  ^  relatly9ly  om(lll  „b„  of  elenenta 

it  is  convenient  to  look  upon  the  vectorial  deteroining  para 
meter  of  a  system  »(t)  an  a  function  of  J***”1" 

ing  parameters  of  elements  H^t),...,  Hbsct)-  At  a 
moment  of  time  t,  the  random  vaotor  Hj  in  a  function  of 
Tandem  vector.^  v,  .  each  of  which  Ln  several  components. 
The  linearisation  of  the  dependence 

•  WH  -  <Pr(HU{,  ....  Hlo{;  H2lt,  •  •  •  >  H«w*»  •  *  -  Hnli’  •  *  M  Hn!t  (3‘i) 


by  expansion  in  a  Taylor  series  and  the  elimination  of 
mi  terms  of  higher  than  the  first  order,  find  th® 
approximate  values  for  the  mathematical  expectations  of 
vector  V/.  components  by  means  of  the  formula 

*  t  tr% 


<pf  («r 


m*iU» 


J  <3.2* 


In  formulas  (3,1)  and  <3, 2)  z  is  the  num 

her  of  components  of  each  of  the  random  vectors  Hli»^|» 

H..}  the  components  of  the  random  vector  W±  are  indica 

#  *  •  f  T>  5  9 

ted  by  the  i rtdex  r  =  t  e »  •  t  ®* 

With  uncorrelated  vectorial  random  processes 
H,(t)  {%  -  1,2, •».,»>  having  uncorrelated  components, 

tfife  dispersions  of  the  random  vector  V±  components  will  be 


«  v  ,  , 
>1  2  (  — 


01ihi  /vnt,M 


in  formula  (3.3)  is  the  dispersion  of  the 

0,-th  component  of  the  random  vector  HiL  i  *L  i*  the  number 
of  components  of  each  of  the  random  vectors  H£. .  Mion 
H  . .H  are  scalar  random  quatitiss,  formula  C3^3J 

11  *  *  *  *  R1 


li *  *  *  *  *  ni 
takes  the  form 


V  {  s%\ 
dtp?;  ■  '  2j  \  (iti.j  ! 


Y  ~)  e«H* 

^  V  <«!»  >‘f„i 
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Having  calculated  the  values  and  Sr  for 

.  A  +  it  is  possible  to  us® 

two  s<mc»ts  °£,t*”a,!ithe  lolt^eotlon  to  fled  the  seen 

the  fowules  Blven  in  t  durmcc  which  the  element 

JTSrSS^SS-^  -  interval  <tt  t  t> 

Assuming  S&’STES-  £  *« 

vectorial  linear  random  function 
mean-square  deviations 

a«V(i+li))  “  COI1&t 

+-  find  the  entire  distribution  law  for  the 
it  is  possible  to  fi  .  order  itt  a  POUgh  approxima- 

time  the  system  is  in  ”  “  ®10B  0l  the  law  of  distribution 
tion.  A  more  s  flna  the  characteristics  mwrCt) 

f<eV*1t>rof  the  components  of  the  vectorial  linear  random 

’r  »(t)  To  calculate  these  characterietics  by  means 
function  W<t).  To  reference  1,  it  is  necessary  to 

of  the  formulas  found^in  ^  ^  raadom  quantities  Wri  and 

Itnow  the  *®up*o  for  the  latter  may  be  obtained  by  sub- 

«r(i+l)'  The.  thG  correlative  moment  formula  values  of 
stituting  into  th  iaed  by  expansion  in  a  Taylor 

s  er  i  es^Af  t  er^  t  he  *  t  rans  f  orraat  i  one  we  have 

’l  "1  /  6>r  \  (  d(?r  \  f.  (3.5) 

A«V iu'r(i+l)  Zi  Zi  V  >  m  V  '  ni,,  t  > 

.  .  {»!  v— -1  ■CvJi  .  .  ,vr” 

_  .  s>  v  K  v  ..  a  is  the  coupling  moment 

In  formula  k3.o>  K * v 

i  «  for  the  ^-th  component  of  the  determining 
of  the  values  for  the  the  i-th  and  (i*l)"th 

parameter  of  the  t-th  element  *» 

cross— sect ions •  with  essentially 

The  calculation  of  *?1  **  £<1+#  J  #ut  by  methods 
non-linear  dependencea  (3.1?  can  o 

described  ?**®'**°®  of  the  durability  of  a  complex 

The  -  -Avt-iatioal  element  boundaries  with 

system  according  to  .  ■ p  is  ln  many  ways  analogous 

vectorial  determining  with  scalar  element  dotor- 

to  the  COTTl*ll*a  Reference  3).  The  difference  consists 
mining  parameters  (  und<jr  consideration  involves  the 

tLdoa  vector*  2^  related  to  moment  of  time  t4  with  components 

ZvS-//vi  —  G.i  (3'6) 


for  \ 

(--) 

V 

rn 

v.on.i  > 

In  formula  (3,5?  h 


8 


In  formula  (3.6)  Hv;  In  the  value  of  the  ^-th 
coxae  one  at  of  tb®  determining  parameter  of  the  element 
at  tin®  U,  and  fl,v-  i®  the  value  of  the  same  component 
of  tb©  vectors al  statistical  bound  of  this  element.  With 
nornal  distribution  laws  and  uncorrelated  random  vector 
components  H4  and ^ .  the  mean  probability  density  of  the 
time  during  which  the  element  is  in  working  order  over  the 
period  (ti ,t.  +i)»  in  accordance  with  formulas  (2,10)  and 
(3»6)«  will  b® 

(3.7) 


Am  'f  U /  is  a  Gaussian  function,  and 


m.n  ,  —  m,.,.  , 
'Jvi 


In  formula®  (3.7)  and  (3.8),  and  are  tbe 

mathematical  expectations  of  the  values  of  the  v -»th 
components  of  the  determining  parameter  and  the  statistical 
bound  at  time  *i?CT^.  and  v;  ar®  the  m.am-square  deviations 

of  these  quantities  from  their  mathematical  expectation®. 

In  order  to  find  the  distribution  law  f(t)  for  the 
time  during  whiolx  the  ©lament  is  in  working  order  j,  it  is 
necessary  to  determine  from  formulas  (1,3)-(1*7>  the  aussori- 
cal  characteristics  which  determine  the  dependences  mY,  (t>, 
a  .<t>,%<*>#  andX>Ct).  Then  formulas  (3.7)  and  (3.8)  are 

used  to  find  the  tjt  values  for  each  interval  which  are  smoothed 

out  by  a  continuous  curve. 

Taking  Into  consideration  correlations  between  the 
components  of  random  vectors  B±  and  ^  considerably  compli¬ 
cate®  th®  process  of  calculating  systems  durability  on  the 
basis  of  statistical  element  bounds  with  vectorial  deter¬ 
mining  parameter®. 
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